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Abstract 

In this paper the first non-linear force-free Vlasov-Maxwell equilibrium is presented. One com- 
ponent of the equilibrium magnetic field has the same spatial structure as the Harris sheet, but 
whereas the Harris sheet is kept in force balance by pressure gradients, in the force- free solution 
presented here force balance is maintained by magnetic shear. Magnetic pressure, plasma pressure 
and plasma density are constant. The method used to find the equilibrium is based on the analogy 
of the one-dimensional Vlasov-Maxwell equilibrium problem to the motion of a pseudo-particle 
in a two-dimensional conservative potential. This potential is equivalent to one of the diagonal 
components of the plasma pressure tensor. After finding the appropriate functional form for this 
pressure tensor component, the corresponding distribution functions can be found using a Fourier 
transform method. The force-free solution can be generalized to a complete family of equilibria 
that describe the transition between the purely pressure-balanced Harris sheet to the force-free 
Harris sheet. 

PACS numbers: 52.20.-j, 52.25.Xz, 52.55.-s, 52.65.Ff 
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Force-free magnetic fields, i.e. magnetic fields satisfying 



(V X B) X B = 0, (1) 
V • B = 0, (2) 

are important for modelling \ow-(3 plasmas in laboratory, space and astrophysical 
applications[l]. Equation ([T]) implies that V x B (basically the electric current density) 
has to be aligned with the magnetic field, i.e. V x B = aB. The scalar function a is 
constant along magnetic field lines due to Eq. ([2]), but can vary from field line to field line. 
If a does not vary from field line to field line, but is globally constant we get the case of 
linear force- free fields (sometimes also called constant-a fields). All other force- free fields 
are called non-linear force-free fields. 

Using magnetohydrodynamics (MHD) many useful linear and non-linear force-free mag- 
netic fields can be found analytically, especially if translational or rotational symmetry of 
the solutions is assumed (see e.g. [ll, y|)- This is completely different if one considers 
collisionless Vlasov-Maxwell (VM) equilibria (see e.g. the discussion in J2|]i. So far, only 
one-dimensional linear force-free VM equilibria have been found0, 0, S, E], and, to the 
best of our knowledge, no non-linear force-free VM equilibria are known. 

One-dimensional (ID) VM equilibria are frequently used as a starting point for stud- 
ies of waves and instabilities in collisionless plasmas. One of the most commonly used 
ID VM equilibria is the Harris sheet 8|, with B{z) = Bota.nh.{z/L)ex and j(z) = 
Bq/^hqL) cosh^'^{z/ L)ey, so the current density is perpendicular to the magnetic field. The 
force balance is maintained by a pressure gradient. Often, a constant magnetic field in the 
y-direction (guide field) is added, which, if sufficiently strong, is used to mimic a force-free 
field. It is clear that through introducing a guide field the current density is partially field- 
aligned, but the strength of the guide field is completely decoupled from the strength of the 
current density. In force-free fields a stronger current density would lead to a stronger shear 
of the magnetic field as the two are closely coupled. Furthermore, a constant magnetic 
field will not add any free energy to the system, whereas one expects an increase in free 
energy if the magnetic shear in a force-free field is increased. As a final point we mention 
that force-free equilibria will have constant density and pressure, whereas the Harris sheet 
plus guide field has the same pressure and density gradients as the Harris sheet itself. This 
may be an important difference in studies of, for example, magnetic reconnection (see e.g. 



0, Q, 
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13|). Some investigations of the stability and dynamics of the known hnear 



force-free ID VM equihbria have been undertaken s], Q, ll5| , but it is to be expected that 
non-hnear force free equihbria will have new and interesting properties. 

Generally VM equilibria can only be found easily for cases with spatial symmetries, and to 
obtain analytical force-free solutions one has to investigate situations with invariance along 
two coordinate directions. In this paper we consider the case of translationally invariant 
VM equilibria depending only on one spatial coordinate, here taken to be z. 

We assume that the magnetic field has components and By. The magnetic field 
components are written in terms of a vector potential A = (A^;, Ay, 0) where B^ = —dAy/dz 
and By = dA^/dz. 

We assume a plasma consisting of two particle species of equal, but opposite charge (elec- 
trons and ions/protons). Due to the symmetries of the system the three obvious constants 
of motion for each particle species are the Hamiltonian or particle energy for each species s, 
Hs = \''^s{v1 + Vy + v1)+qs(t), the canonical momentum in the x-direction, p^s = ^s^x + ^ls^xi 
and the canonical momentum in the ^/-direction, pys = nisVy + qsAy. Here (p is the electric po- 
tential and rris and are the mass and charge of each particle species. All positive functions 
fs satisfying the appropriate conditions for existence of the velocity moments and depending 
only on the constants of motion, fg = fs{Hs,Pxs,Pys) are solutions of the steady-state Vlasov 
equation. 

One can show 16|, ll7|] that for a quasi-neutral plasma. Ampere's law can be written as 

d^A^ _ dP,, 

~d^ - ~^'^x' ^ ^ 

d^Ay dPzz 

-d^ = 

where Pzz{Ax, Ay) is the zz-component of the plasma pressure tensor, defined by 

/oo 
msV^Jsd\. (5) 

Equations ([3]) and (jl]) can be immediately integrated once to give the force balance condition 
across the sheet as 

1- Pzz = Pt = constant. (6) 

2/io 

Due to Eq. ([T]) a force-free equilibrium satisfies the conditions B^ = constant and Pzz = 
constant separately. 
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FIG. 1: The magnetic field, current density and pressure profiles as functions of zjL for the Harris 
sheet (left panel), the force- free Harris sheet (right panel) and intermediate case (middle panel). 

The ID VM equilibrium equations ([3]) and (jlj) are equivalent to the equations of motion of 



a (pseudo-) particle in a conservative 2D pseudo-potential Pzz{.-^x-, -^y)^ ll6|. The position 
of the pseudo-particle is given by A^^ Ay with the pseudo-time given by z. The energy 
(Hamiltonian) of this pseudo-particle is given by the total pressure defined in Eq. ([6]) 
(modulo a factor /io) E = [{dAx/dz)"^ + {dAy/dzY]/2 + fioP^ziA^, Ay). One can show that 
a force-free VM solution corresponds to a pseudo-particle trajectory that is identical to a 
contour of the pseudo-potential [l7|]. This is easily possible for attractive central potentials 
which have circular contours and also allow circular pseudo-particle orbits. These circular 



orbits correspond to the known linear force- free solutions j^. 



7|, which, as far as we are 



aware, are the only known force-free VM solutions. For finding nonlinear force-free solutions 
we obviously need to find a pseudo-potential {P^ziA^, Ay)) which is not a central potential, 
but still allows a solution to Eqs. ([3]) and (j4]) that is identical with an equipotential line. 

Channellj^ showed how, by making a number of sensible assumptions, a transform 
method can be used to determine a class of distribution functions for a known Pzz{A^, Ay). 
Mynick et al.|l6| generalized this method and used it to determine the distribution functions 
numerically. In this paper we will first determine a function Pzz{Ax, Ay) for the force-free 
Harris sheet and then use Channell's method to find the corresponding distribution func- 
tions. 

The magnetic field of the force-free Harris sheet solution is given by (see also Fig. [1]) 

= Botanh{z/L), (7) 
By = Bocosh~\z/L), (8) 

with Bq the constant amplitude of the field and L the sheet half width. Obviously we 



have Bl + By = Bq. One can easily see that jx = Bq/ (fioL) tanh{z/ L)/ cosh(z/L), jy = 
Bq / {^qL) co^hr'^ {z / L) , giving a{z) = [L cosh(2;/L)]"^. 

The x-component of this magnetic field is identical to the Harris sheet B^, but in this 
case the force balance is maintained by the magnetic shear component By instead of the 
plasma pressure. The vector potential for the force-free Harris sheet field is found to be 
given by 

^xjfh = 2BoL aictan {exp{z / L)) , (9) 
Ayjfh = -BoLln (^cosh (^|)) , (10) 

in a convenient gauge. 

In order to make analytical progress we assume that Pzz has the form Pzz{Ax, Ay) = 
Pi{Ax) + P2{Ay). The physical meaning of this assumption is that for each particle species 
there are two different particle populations that carry the components of the current density 
in the x- and the ^/-directions. Eqs. ([3]) and @ give the conditions 

2 

I +2/ioPi(A,) = 2/ioPoi, (11) 

2 

] + 2fioP2{Ay) = 2/ioPo2, (12) 



dA^ 
dz 

'dA. ' 



dz 



where Pqi and P02 are constants. Equations (fTT!) and (fT2l) will be used to find the appropriate 

Pzz{Ax, Ay). 

We substitute A^jfh and Ayjf^ into the first terms of Eqs. fllip and fll2p and then use 
that exp{z/L) = tan(y4^./2Po-^) and cosh(2;/L) = exp{—Ay/BQL) to obtain 

Pg 1 / 2Aj. \ / 2 Ay 



P. 



2/io 



- cos + exp , 

2 \BoLj ^ \BoL 



+ Po3. (13) 



A surface plot of Pzz{Ax, Ay) is shown in Fig. [21 Above the surface plot the trajectory 
representing the force-free Harris sheet solution in the Ax-Ay-p\a.ne is shown. By construction 
it is identical to a contour of PzziA^, Ay). 

We use Channell'slsl Fourier transform method to solve the integral equation ([5]) for the 
distribution functions fg. The method is based on the assumptions that (a) the distribu- 
tion functions have the form fs{Hs, Pxs,Pys) = fosQW{-PsHs)gs{Pxs,Pys) and that (b) the 
quasineutral electric potential (j)qn vanishes. The validity of the second assumption can be 
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FIG. 2: A surface plot of the pressure function Pzz{Ax, Ay) for the force-free Harris sheet. The 
force- free Harris sheet solution is shown as a pseudo-particle trajectory at the top of the plot. It 
is identical with a contour of Pzz- 

easily verified a posteriori and only requires the correct choice of parameters. Applying the 
method we find that the required distribution functions are of the form 



fs 



'^th,s 



exp{-f3sHs) [exp (PsUysPys) 



+as cos iPsUxsPxs) + bs 



(14) 



where Vth,s = ("^s/^s)"^''^ is the thermal velocity of particle species s and u^s, Uys, as and bg 
are constants with < < 6^- We have here reverted to the usual microscopic notation 
for the distribution functions. We will make the connection to the notation used previously 
by calculating P^z directly from the distribution function and then comparing the result 
with the expression (fT3l) . This is useful to relate the macroscopic quantities Bq and L to 
the microscopic parameters of the distribution function. The first part of this distribution 
function is identical with the Harris sheet Sj distribution function, whereas the second part 
corresponds to a different particle population which carries the current in the x-direction 
and is responsible for the shear field By{z). When calculating Pzz from the distribution 
function one finds that it has the general structure Pzz = J2s ^wi~QsPs(p)Nsi^x, ^y), 
where 



+ as exp{-Psms{uls + cos^PsUxsgsA^) 
+ bsexp{-Psmsuls/2)]. 



The quasi-neutrality condition leads to 0g„ = [e(/3e + ^^{Ni / Ne) ■ The condition of 

vanishing quasi-neutral electric potential implies that Ni{Ax,Ay) = Ne{A;j:, Ay), which is 
true if 



aeexp{-/3eme{ul^ + Wye)/2) = 

tti exp{-Pimi{uli + uli)/2) = a, 
6eexp(-/3emeMyg/2) = biexp{-(3imiulJ2) = b, 

Supposing that f3e and /5j are given we have ten other parameters needing to satisfy only five 
equations, which is always possible. This provides the necessary a 2)ostenon justification for 



using Channell's method. Using the notation often used for the Harris sheet (e.g. 18|) P^z 
becomes 

Pzz = + no[exp{-e(3eUyeAy) 

+ a cos(e/9eM^e^z) + b] . (15) 

Comparison with Eq. (JT3l) shows that for the force-free Harris sheet the connection 
between the microscopic notation and the original notation is given by 

-^i + D no, (16) 



2/io \Pe A 

^/ioe2noMye/5e(/?e + A) 



« = (18) 
b = 2fioPo3/Bl (19) 

Eq. f|T7|) is especially important as it provides a relation between the length scale L and 
the parameters of the distribution function. This, for example, directly links a{z) = 
[L cosh(2;/L)]~^ derived from the general form of the magnetic field and current density 
to the microscopic parameters of the distribution function. 

It is straightforward to see that for different parameter values the distribution function 
( fT^ gives the complete family of equilibria describing the transition between the Harris 



sheet and the force-free Harris sheet. The intermediate cases have, written as functions of 



= S,otanh(2/L), (20) 
By = Byocosh-\z/L), (21) 
P,, = PoCosh-^(^/L)+Poo, (22) 



where Pq + ByQ/2fiQ = P^Q/2/io. Taking the hmit ByQ gives the Harris sheet |8| and 
taking the hmit Pq ^ gives the force-free Harris sheet. 

An appropriate P^ziA^, Ay) can be determined in the same way as for the force- free Harris 
sheet and takes the form 

Bin / 2A„ \ 1 B^n ( 2A^ \ , , 

P,^ = ^exp ^ + ^cos ^ +Po3. (23) 



In this case a comparison between Eqs. fl23l) and f|T5l) shows that Eqs. f|T6l) . f|T7|l and ([19 
do not change apart from Bq becoming P^-o, but that we now also have 



-1/2 



As shown by Harrison and Neukirch|17l] one can deduce from one Pzz{Ax, Ay) allowing a 
force-free VM solution an infinite number of other functions Pzz{Ax, Ay) allowing the same 
force-free solution by using positive definite functions of the known Pzz{Ax, Ay). We mention 
in particular the possibility of using an exponential function of the Pzz presented here, as 
it would give rise to a product form for Pzz instead of a sum. The distribution functions 
would also consist of products of functions of p^s and Pys instead of a sum. It is, however, 
unclear whether the method used in this paper would still allow for an analytical calculation 
of these distribution functions. 

This new family of VM equilibria will generate new possibilities for studies of linear 
and nonlinear instabilities of force-free current sheets. The stability of the VM equilibria 
presented here have yet to be investigated. We point out that the p^^s-dependent part of the 
distribution function may have multiple peaks in the f^^-direction and we suspect that this 
will give rise to instabilities. We also remark that although the Bx{z) and jy(2;)-profiles are 
identical to the Harris sheet, jx{z) is antisymmetric with respect to z = 0. This is closely 
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linked to the fact that in the Harris sheet solution the spatial structure of the current 
density is determined by the density structure with the average velocity of each particle 
species being constant, whereas in the force- free solution presented here the particle density 
is constant and the spatial structure of the current density is determined by the spatial 
structure of the average velocity. Further investigations will be needed to clarify exactly 
what the implications are for the stability of the new solution, but on the basis of the 
physical differences just mentioned one would expect the stability properties of the force- 
free solution to differ considerably from those of the Harris sheet. Apart from studying the 
stability properties of the solution class presented here, it will be also be very interesting to 
investigate whether the general method employed here can be used to find other non-linear 
force-free solutions. 
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